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1.  INTRODUCTION 


In  the  National  Airspace  System  (NAS),  the  corrected  slant  range  R(S,H)  of  a  target 
relative  to  a  radar  is  determined  from  the  measurement  S  of  slant  range  and  the 
reported  altitude  H.  The  computation  of  R(S,H)  involves  the  calculations  of  a 
square  root  function  using  an  approximation  technique.  The  present  approximation 
was  developed  at  a  time  when  the  only  use  made  of  radar  surveillance  data  was  for 
display  purposes  involving  movement  of  the  data  block  on  the  planned  view  display 
The  computational  accuracy  requirements  for  this  purpose  were  of  little  signifi¬ 
cance.  Evolution  of  the  air  traffic  control  system  towards  a  greater  level  of 
automation  (in  providing  advanced  air  traffic  control  functions  such  as  conflict 
alert)  now  places  substantially  increased  requirements  on  the  accuracy  of  the 
functions  which  support  the  radar  cracking  algorithm.  Unfortunately,  the  need  for 
increased  computational  accuracy  in  support  of  the  advanced  automation  features  has 
been  recognized  only  recently.  It  is  the  objective  of  this  report  to  provide  an 
approximation  technique  with  a  specific  predetermined  error  boundary  which  will 
guarantee  that  the  computational  errors  resulting  from  the  function  approximation 
will  have  absolutely  no  measurable  impact  whatsoever  on  any  aspect  of  system  per¬ 
formance.  It  is  important  for  whatever  technique  is  employed  that  the  algorithm 
be  adaptable  to  the  accuracy  of  the  present  radars  as  well  as  those  in  the 
immediate  future. 

If  the  radar  site  is  located  at  an  altitude  HR  above  mean  sea  level,  then 
(reference  1) 


R(S,H) 


=  (s2  -[h-hr]2)1/2 


(1) 


Due  to  the  necessity  of  performing  such  computations  in  real  time  for  numerous 
targets,  an  approximation  of  the  square  root  function  is  employed  in  NAS  that  does 
not  severely  tax  available  computational  resources.  A  previous  investigation  of 
this  approximation  (reference  2)  demonstrated  that  errors  as  large  as  0.3  nautical 
miles  (nmi)  can  be  introduced  in  the  corrected  slant  range.  Also,  it  is  shown 
there  that  an  adjustment  of  the  parameters  used  in  the  current  approximation 
algorithm  would  eliminate  errors  in  excess  of  0.125  nmi,  the  present  quantization 
of  slant  range  measurements.  However,  this  parameter  change  will  not  reduce  the 
computation  time  for  slant  range  correction.  Also,  because  of  the  structure  of 
the  present  FAA  algorithm  it  is  not  easily  adaptable  to  an  increased  approximation 
accuracy  specification.  This  report  deals  with  an  alternative  to  the  approximation 
in  current  use.  It  can  be  adapted  to  meet  any  reasonable  specification  of  the 
maximum  allowable  approximation  error. 

Briefly  stated,  our  problem  is  as  follows.  The  altitude  HR  of  the  radar  site  is 
a  known  constant,  S  is  observed  at  the  radar  site,  and  H  is  observed  at  the  target 
and  then  transmitted  to  the  ground  in  response  to  an  appropriate  interrogation 
signal  from  the  radar.  In  addition,  there  is  a  well-defined  region  of  admissible 
values  for  S  and  H,  that  is  determined  by  the  measurement  capabilities  of  the  radar 
and  the  altitudes  of  targets  of  interest.  The  objective  is  to  find  a  simple 
estimate  R  (S,H)  of  R(S,H)  such  that  the  absolute  value  of  the  difference  between 
the  two  does  not  exceed  a  predetermined  error  bound  c  over  the  admissible  region. 
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The  region  of  admissible  values  for  Che  measurement  vector  (S,H)  is  defined  in 
section  2.  Basic  elements  of  the  algorithm  and  machine  implementations  of  the 
algorithm  are  described  in  section  3. 

Section  4  describes  how  to  choose  the  parameters  of  the  algorithm  so  that  the 
approximation  error  will  never  exceed  the  predetermined  error  bound  e  over  the 
admissible  values  of  S  and  H.  The  algorithm  is  illustrated  by  a  numerical  example 
in  section  5,  and  comparisons  are  made  with  the  estimation  procedure  currently 
employed  in  NAS.  Section  6  contains  possible  modifications  of  Che  algorithm  and 
section  7  contains  conclusions. 


2.  CONSTRAINTS  ON  SLANT  RANGE  AND  ALTITUDE 


There  exist  combinations  of  reported  altitude  and  slant  range  measurement  that  are 
unacceptable  in  the  determination  of  corrected  slant  range.  For  example,  admis¬ 
sible  values  of  S  are  limited  by  the  maximum  effective  range  S^  of  the  radar  as 

well  as  interference  phenomena  encountered  at  distances  less  than  some  minimum 
range  Sm.  In  addition,  due  to  limitations  on  altitudes  of  targets  of  interest,  H 
can  be  assumed  to  be  bounded  above  by  some  constant  K.  Moreover,  the  ratio  of 
|H-Hr|  to  S  is  constrained  by  the  so-called  cone  of  silence  of  the  radar  to  be 

less  than  another  constant  J.  In  other  words,  S  and  K  can  be  viewed  as  being 

restricted  under  practical  operating  conditions  by  the  following  inequalities: 


0<sm  <  s  <  SM 

(2) 

h-Hr|  <k 

(3) 

H-Hr|  /S<J<1 

(4) 

In  a  typical  NAS  air  route  traffic  control  center,  radar  site  altitudes  are  less 
than  1.646  nmi  (10,000  feet),  SM  does  not  exceed  200  nmi ,  Sm  is  at  least  2  nmi ,  J 
is  sin  70°,  and  controlled  traffic  consists  mainly  of  aircraft  at  altitudes  less 
than  9.875  nmi  (60,000  feet).  Thus, 

K/Sm  <J  (5) 

in  the  case  of  NAS,  and  this  inequality  will  be  assumed  to  be  satisljd  throughout 
the  remaining  discussion. 


3.  ALGORITHMIC  STRUCTURE 


The  reported  altitude  H  and  the  measured  slant  range  S  are  hereafter  assumed  to  be 
constrained  by  relations  (2)  -  (4)  where  HR,  SB,  S^,  K,  and  J  are  known  constants. 
Subject  to  these  constraints,  our  objective  is  to  find  a  simple  approximation 
R,(S,H)  of  R(S,H)  such  that 
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j R  ( S , H)  -  R(S,H)|  <  C 


(6) 


where  e  is  some  prescribed  error  bound,  say  0.125  nmi .  In  what  follows,  it  will  be 
shown  how  this  desideratum  can  be  accomplished  by  an  algorithm  based  upon  a  parti¬ 
tion  of  the  interval  between  0  and  J2  into  a  finite  number  n  of  subintervals  by 
boundary  points 


o=y.<y  »< • • -<yn  =  J2 


(7) 


The  algorithm  is  of  the  following  form:  if 


yk-i  £<Ih_hrI  /s)2<yk 


(8) 


where  It  is  any  one  of  the  integers  1  through  n,  then  set 
Ro(S,H)=S  j  Ar  +  \  (|h-hr|  /S)2 


(9) 


The  problem  is  to  find  n,  the  n-1  interior  boundary  points  and  the  2n 
coefficients  Aj,...,An,  As  will  be  seen,  these  depend  upon  the 
error  bound  e.  Moreover,  it  will  be  shown  how  these  parameters  can  be  selected  so 
that  Ro  (S,H)  is  never  less  than  R(S,H).  Thus,  if  instead  of  R^  (S,H),  the 
corrected  slant  range  is  approximated  by  Ro(S,H)-e/2,  then  the  approximation  error 
will  be  e/2  rather  than  e-  In  fact,  it  will  be  shown  that  if  the  latter  approxi¬ 
mation  is  used  as  an  initial  estimate,  then  after  p  applications  of  the  Newton- 
Raphson  iteration  (reference  3),  the  corrected  slant  range  can  be  found  to  within 
an  error  of  (e/2)2  • 


Before  proceeding  with  the  mechanics  of  parameter  selection  for  the  approximat ion 
algorithm,  its  implementation  will  be  discussed.  First  of  all,  we  point  out  that 
the  ratio  of  |H-Hr|  to  S  is  determined  in  NAS  for  the  purpose  of  controlling  bias 
errors  in  radar  measurements.  Thus,  only  the  square  of  the  ratio  is  unique  to  the 
algorithm.  However,  we  will  show  that  it  is  possible  to  eliminate  the  squaring 
operation  at  the  expense  of  memory.  On  the  other  hand,  as  will  be  seen  shortly, 
the  difference  S2  -(H-Hr)2  is  required  by  the  Newton-Raphson  iteration. 
Hence,  if  the  approximation  algorithm  is  to  be  used  to  supply  that  iterative 
technique  with  an  initial  estimate  of  the  corrected  slant  range,  then  there  is 
little  reason  to  oicrifice  memory  to  eliminate  an  operation  that  can  be  efficiently 
employed  to  determine  a  necessary  element  involved  in  the  calculation  of  the  final 
resul t . 


The  second  implementation  item  to  consider  is  the  obvious  need  for  some  procedure 
to  identify  the  particular  integer  k  for  which  the  relation  (8)  is  satisfied.  This 
can  be  accomplished  by  a  linear  search.  Another  possibility  is  to  subdivide  the 
continuum  from  0  to  1  into  intervals  of  identical  length,  equal  to  some  negative 
power  of  2,  and  employ  a  hashing  technique  in  which  higher  order  bits  of  a  fixed 
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point  representation  of  ( J  K—  H r [ / S ) 2  are  used  to  make  the  appropriate 
identification.  For  example,  suppose  m  is  any  integer  for  which  2~m  does  not 

exceed  J2  f  and  let  N(ra)  represent  the  smallest  integer  j  for  which  j2  111  exceeds  or 
is  equal  to  J2.  Then  the  sets  defined  by 


U 


J 


{t : ( j-1 )2~™  <t  <  j2~m} 
{t:[N(m)-lJ  2~m<t<J2} 


if  j  =  1 . N(m)-1 

if  j=N(m)  (10) 


partition  the  continuum  from  0  to  J2  into  N(m)  disjoint  intervals.  Suppose  that 
(JlH-Hgp/S)2  is  a  member  of  U j  •  Then  the  m  highest  order  bits  in  a  fixed  point 
representation  of  that  ratio  constitute  the  binary  representation  of  j-1,  and  this, 
in  turn,  is  an  automatic  signal  to  the  effect  that  the  ratio  is  indeed  a  member  of 
Ui.  Now  choose  m  to  be  sufficiently  large  that  each  interval  contains  no  more 
than  one  of  the  boundary  points  y»,-».,yn;  e.g.,  choose  m  large  enough  to  satisfy 
the  relation 


2_m  <min  (yk  -y^) 

1<  k<_  n 

In  the  event  that  Uj  does  not  contain  any  one  of  the  interior  boundary  points 
yi,...,yn_i,  our  binary  signal  is  tantamount  to  an  automatic  identification  of 
the  integer  k  for  which  (8)  is  satisfied.  Otherwise,  a  collision  occurs  in  the 
sense  that  Uj  contains  both  the  ratio  and  one  of  the  interior  boundary  points.  In 
this  situation,  we  are  confronted  with  one  of  the  possibilities;  i.e.,  the 
offending  boundary  point  is  identical  to  the  smallest  element  (j-l)2-m  in  the 
interval,  or  else  it  exceeds  this  number.  If  it  exceeds  the  minimum  element,  then 
it  must  be  compared  with  the  ratio  in  order  to  make  the  correct  identification  of 
the  integer  k  for  which  (8)  is  satisfied.  If  it  is  equal  to  the  minimum  element, 
then  the  relative  positions  of  the  ratio  and  the  boundary  point  are  known,  and 
there  is  no  need  for  a  comparison.  Needless  to  say,  the  likelihood  of  a  collision 
decreases  as  m  increases.  In  this  way,  computational  speed  can  be  increased  at  the 
expense  of  memory. 

The  hashing  technique  that  we  have  described  will  never  require  a  comparison  if  the 
minimum  element  of  each  interval,  but  the  first  is  an  interior  boundary  point; 
i.e.,  n*N(m)  and  yj  *  j2-m  for  all  j*l,...,n-l.  As  will  be  shown  later,  subject 
to  the  condition  that  the  minimum  elements  of  the  Uj's  meet  this  requirement,  there 
is  a  way  to  pick  the  smallest  number  m  of  intervals  such  that  the  approximation 
error  satisfies  (6).  In  other  words,  without  sacrificing  an  undue  amount  of 
memory,  we  will  show  how  to  subdivide  the  continuum  from  0  to  J2  in  such  a  way 
that  collisions  do  not  require  comparisons. 


4.  PARAMETER  SELECTION 


Turning  now  to  the  problem  of  parameter  selection,  we  find  it  convenient  to  work 
with  the  independent  variables  S  and 
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x  =  ( I H-Hr I /s ) 2 


(12) 


rather  than  S  and  H.  Obviously,  the  corrected  slant  range  (1)  can  be  expressed  as 
a  function  of  these  variables,  namely, 


R  (S,  x)  =  S  (l-x)1/2 


(13) 


Likewise, _ the  estimate  (9)  provided  by  the  algorithm  can  be  represented  as  a 
function  R„  (S,x)  of  the  same  variables.  Also,  the  constraints  (2)  -  (4)  can  be 
expressed  in  terms  of  equivalent  restrictions  on  S  and  x.  In  particular,  letting 
xA//  represent  the  positive  square  root  of  x,  (3)  implies 


S  =  |h-Hr|x-1/2  <Kx-1/2 


(14) 


On  the  other  hand,  the  right  side  of  (14)  exceeds  S^  only  if  x  falls  below  (K/S^)2. 
From  this  we  conclude  that  the  constraints  (2)  -  (4)  are  equivalent  to  the  relations 


0  ^  x  <  J2 

(15) 

<  Sm  <S  <SM 

if  0  <x  <  (K/Sm)2 

(16) 

<  Sm  <S  <K*~1/2 

if  (K/Sm)2  <x  <  J2 

(17) 

It  now  remains  to  show  how  one  determines  the  estimate  R  (S,x)  of  R(S,x)  so  that 

o 

the  absolute  difference  between  the  two  does  not  exceed  e  for  all  vectors  (S,x) 
satisfying  relations  (15)  -  (17). 


Suppose  0  is  a  nonnegative  number  less  than  1.  Then  the  tangent  line  to  the  graph 
of  (l-x)l/2  at  x= 0  can  be  represented  by  the  function 


where 

and 


a(x,0)  =  A(0)  -  B(0)  x 
A(0)=(l-0/2)(l-0)"1/2 
B(0)=(1-0)~1/2  12 


(18) 

(19) 

(20) 


Since  the  second  derivative  of  (l-x)2/^2  is  negative,  it  follows  that 

e(x, 8)*a(x, 0)  -  (1-x)1^2  >  0  (21) 


for  all  x  in  the  continuum  from  0  to  1.  Hence,  if  corrected  slant  range  (13)  is 
approximated  by  Sa(x,0),  then,  after  multiplying  (21)  through  by  S, 
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Se(x.B)  =  S  a  (x,fl)  -  R(S,x) 


(22) 


is  the  error,  and  it  is  always  nonnegat i ve. 
(17)  that 

0  <Se(x,B)  j<  E  (x,e) 


Moreover,  it  follows  from  (16)  and 

(23) 


where 


E(x,6) 


{^Me 

Ke(x , B)x'1/2 


if  0  <x<(K/Sm)2 
if  (K/Sm)2  _<  x  <J2 


(24) 


In  fact,  it  is  clear  that  E(x,8)  is  the  maximum  error  that  can  be  incurred  through 
use  of  the  approximation  Sa(x,$)  so  long  as  S  and  x  satisfy  the  constraints 
(15)  -  (17). 

It  can  be  verified  by  differentiation  that  when  E(x,B)  is  viewed  as  a  function  of  x 
alone,  it  decreases  monotonical ly  as  x  increases  from  0  to  B  >  it  vanishes  at  x-  g , 
and  it  increases  monotonical ly  as  x  increases  from  B  to  1.  Thus,  corresponding  to 
any  positive  value  that  one  might  care  to  assign  to  c,  there  must  exist  numbers 
xj(B)  and  X2(0)  such  that 


and 

so  long  as 


0  <X!  UK  B  <x2 

(25) 

E(x,bK  e 

(26) 

XiU)  j<  x  <  X2  (fl) 

(27) 

The  situation  is  illustrated  in  figure  1.  In  other  words,  if  the  restriction  (27) 
is  added  to  the  constraints  (15)  -  (17),  then  S  a(x,B)  approximates  the  corrected 
slant  range  to  within  an  error  e-  Thus,  from  figure  1,  it  appears  that  by  picking 
several  different  B's  between  0  and  1,  we  should  be  able  to  develop  an  algorithm 
along  the  lines  of  (7)  -  (9)  that  will  estimate  the  slant  range  correction  to 
within  an  error  c  for  all  possible  values  of  x  between  0  and  J2. 

The  assertion  at  the  end  of  the  preceding  paragraph  can  be  established  in  a 
rigorous  fashion.  In  particular,  let  us  consider  the  case  where  the  number  n  of 
boundary  points  is  identical  to  N(m) ;  i.e.,  the  smallest  interger  by  which  2  m  can 
be  multiplied  to  yield  a  number  greater  than  or  equal  to  J2,  and  the  interior 
boundary  points  are  given  by 


yj  "  j2"®;  j  *  !•»•••  >n-i 


(28) 


It  is  shown  in  appendix  A  that  for  each  k*l,...,n  there  is  one  and  only  one  number 
ftk  between  0  and  J2  such  that  E(yk_j,  Bk)  and  E(yk,0g)  are  the  same.  Moreover, 
it  is  shown  that  there  exists  a  positive  integer  b<  such  that 
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FIGURE  1.  ERROR  DUE  TO  THE  APPROXIMATION  a(x;0) 


e  (yk 


8  k)  <_  e  for  all  k-1 , .  .  .  ,n 


(29) 


if  10=111°,  and  in  this  case 

E  (x,&kK  e  whenever  yk_j  <_  x  <  yk  (30) 

for  all  k  =  l,...,n.  Now  let  us  choose  the  coefficients  Ap  . .  .  ,  An ,  Bp...,  Bn  of 
the  algorithm  (7)  -  (9)  in  accord  with  the  relations 


Ak  =  A  (0k)  and  Bk  =  B(8k) 


(31) 


where  A(  B )  and  B(B)  are  defined  by  (19)  and  (20).  Then,  when  (  |H-Hk  |  /S)2  is 
substituted  for  x,  the  estimate  R°(S,H)  provided  by  the  algorithm  is  just 

n 

R°  (S , x)=  Sa(x,Bk)  l(x,8k)  (32) 

k=l 

where  a(x,0)  is  given  by  (18)  and 


l(x, 0k)= 


if  Vk-l  ix  <yk 


otherwi se 


(33) 


It  remains  to  show  that  (32)  approximates  the  corrected  slant  range  (13)  to  within 
an  error  e . 

Using  (21)  -  (23)  and  (32),  it  follows  that  the  difference  R#(S,x)  -  R(s,n) 
between  the  estimate  and  the  corrected  slant  range  must  satisfy  the  inequality 


0  <Ro  (S,x)  -  R( s , x)  <E,(x)  (34) 


E.  (x)  “2  E(x,8^)  1(x*6k)  (35) 

k-i 

Moreover,  since  (23)  holds  for  all  S  and  x  satisfying  (15)  -  (17),  the  same  is  true 
for  (34).  Needless  to  say,  if  m*m> ,  then  (30),  (33),  and  (35)  imply  that  Eo(x)  can¬ 
not  exceed  the  error  bound  e  for  all  x  in  the  continuum  from  0  to  J2.  In  other 
words,  when  m*m#  ,  the  estimate  (32)  provided  by  the  algorithm  cannot  be  less  than 
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~  u. 


’v  r  *:• 


the  corrected  slant  range,  nor  can  it  exceed  the  corrected  slant  range  by  an  amount 
greater  than  e  .  Thus,  the  estimate  Ro(S,x)  -  e/ 2  must  approximate  the  corrected 
slant  range  to  within  an  error  e/2. 

By  now  the  reader  should  be  cognizant  of  a  systematic  approach  to  the  design  of  an 
algorithm  for  estimating  corrected  slant  range  subject  to  a  prescribed  limitation 
on  the  estimation  error.  In  particular,  starting  with  m=l ,  one  determines  for  each 
k=l , . . . ,N(m)-l  the  number  Bg  for  which  E( [ k-1 J 2-m(Bg)  and  E(k2-m)  8^)  are  equal  as 
well  as  the  number  8fj(m)  f°r  which  E(  [  N(m)-1  J  2_mi  B^m) )  and  E  ( ,  8(j(m)  ^  are  t*ie 
same.  If  E( [ k-1 J 2_m, 8g)  does  not  exceed  e  for  each  k= 1 ,  .  .  .  ,N (m) ,  then  m  is 
identical  to  mo  .  In  this  case,  one  sets  n=N(m) ,  chooses  the  coefficients 
A1 . A  ,B1,...,Bn  in  accord  with  (31),  and  selects  the  interior  boundary  point  yk 

to  be  k2-m  for  each  k=l,...,n-l.  Otherwise,  m  is  increased  by  1,  and  the  entire 
procedure  is  repeated. 


5 .  EXAMPLE 


We  point  out  that  the  integer  m#  is  dependent  on  the  error  bound  c  .  In  fact, 
as  e  decreases  toward  0,  m  increases  without  bound.  Thus,  the  cost  in  terms  of 
memory  of  implementing  the  algorithm  becomes  prohibitive  as  c  approaches  0. 
However,  this  does  not  mean  that  the  algorithm  has  no  practical  value.  For 
example,  in  current  NAS  operations,  slant  range  measurements  are  quantized  into 
units  of  1/8  nmi ,  and  altitude  is  reported  in  units  of  0.0165  nmi  (100  feet). 
Suppose  we  choose  e  to  be  1/8  nmi.  Then,  under  the  assumption  that  Sm  =  2  nmi, 
SM  =  200  nmi  ,  J  =  sin  70°,  and  K  =  9.875  nmi  in  accord  with  the  operation  of  a 
typical  air  route  traffic  control  center,  it  turns  out  that  =  2  and 


B1  =  0.0630  ,  82  =  0.3705,  B3  =  0.6294,  B4  =  0.8214 

The  coefficients  of  the  algorithm  can  be  determined  from  (19),  (20),  and  (31). 
From  (34),  the  minimum  error  incurred  through  use  of  the  estimate  (32)  provided 
by  the  algorithm  is  0,  and  the  maximum  error  is  E  (x).  As  shown  in  figure  2, 
Eo(x)  never  exceeds  the  prescribed  error  bound  of  1/8  nmi.  Hence,  by  subtracting 
1/16  nmi  from  the  estimate  supplied  by  the  algorithm,  it  is  possible  to  determine 
the  corrected  slant  range  to  within  an  error  of  no  more  than  1/16  nmi. 

A  bimodal  algorithm  is  currently  employed  by  NAS  to  determine  corrected  slant 
range.  In  the  case  where  the  ratio  |H-Hg|/S  is  greater  than  or  equal  to  6/10,  the 
algorithm  computes  and  (H-Hg)^,  and  then,  by  means  of  a  linear  search  procedure, 
it  applies  a  continuous  approximation  consisting  of  six  straight  line  segments  to 
the  difference.  When  the  ratio  is  less  than  6/10,  the  algorithm  approximates  the 
corrected  slant  range  by  the  expression  -  (H-Hg)^  /2S  which  is  equivalent  to  the 
product  of  S  and  the  tangent  to  the  graph  of  (l-x)*'^  at  x  *  0.  Obviously,  both 
modes  of  the  algorithm  require  division  and  squaring  operations.  Moreover,  using 
the  values  assigned  to  Sm,  SM,  J,  and  K  in  our  example,  it  can  be  shown  that 
errors  in  excess  of  1/8  nmi  are  possible  within  the  constraints  (2)  -  (4). 
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FIGURE  2.  ERROR  DUE  TO  THE  APPROXIMATION 


Although 
for  many 
def i ni te 
speed . 


the  bimodal  algorithm  has  met  the  needs  of  the  air  traffic  control  system 
years,  it  is  evident  that  the  aleorithm  of  our  example  offers  some  very 
advantages  in  terms  of  accuracy,  utilization  of  memory,  and  computational 


6.  MODIFICATIONS 


As  already  indicated,  the  amount  of  memory  required  by  machine  realizations  of  the 
approximation  (32)  increases  indefinitely  as  the  error  bound  e  decreases  toward  0. 

To  be  more  specific,  in  our  example,  it  has  been  shown  that  if  e  =  1/8  nmi ,  then  a 
partition  consisting  of  four  intervals  is  sufficient  to  maintain  the  approximation 
error  incurred  through  use  of  the  estimate 


Z  =  R  (S,x)  -  e/2 

Q  O 


(36) 


to  within  e/2  nmi.  However,  if  e  is  reduced  to  1/16  nmi,  then  the  reader  can 
verify  that  a  partition  of  similar  structure  must  have  at  least  8  intervals  in 
order  to  guarantee  that  the  absolute  difference  between  (36)  and  the  corrected 
slant  range  (13)  will  not  exceed  c/2  nmi.  On  the  other  hand,  for  a  given  e,  we  can 
always  use  (36)  as  an  initial  estimate  of  the  corrected  slant  range.  This,  in 
turn,  can  be  used  to  find  a  better  estimate  by  means  of  some  iterative  technique. 

For  instance,  let  us  use  a  to  denote  the  corrected  slant  range  R(S,x).  This  is 
just  the  positive  solution  to  the  equations  obtained  by  setting  w  of 


w 


(37) 


equal  to  0.  According  to  the  Newton-Raphson  iteration,  if  zk  is  the  kth  (K>_0) 
estimate  of  the  solution,  then  the  (k+I)th  estimate  is  just  the  intersection  of 
the  z-axis  with  the  tangent  to  the  graph  of  (37)  at  zmz^,  i.e., 


zk+l  “  (zk  +  <*2/zk)  ,2 


(38) 


In  appendix  B  it  is  shown  that  if  the  minimum  range  Sm  of  the  radar  is  at  least 
2  nmi  and  J  is  sin  70’,  as  is  the  case  in  NAS,  then,  after  p(p>l)  applications  of 
(38)  starting  with  the  initial  estimate  zo  the  error  satisfies  the  relation 

|zp-«|<  (s/2)2P  (39) 

provided  that  the  error  bound  c  does  not  exceed  0.368  nmi.  In  our  numerical 
example,  e  is  1/8  nmi,  and  so  the  estimate  z%  obtained  by  way  of  the  algorithm 
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(7)  -  (9)  must  be  within  1/16  nmi  of  a-  Thus,  after  one  application  of  the 
Newton-Raphson  iteration,  we  are  assured  that  the  absolute  value  of  the  difference 
between  zj  and  the  corrected  slant  range  will  be  no  more  than  1/256  nmi. 

Our  development  of  the  algorithm  (7)  -  (9)  is  based  upon  the  idea  that  one  can  use 
a  tangent  line  to  approximate  a  function  to  within  a  prescribed  error  bound  e  over 
a  small  subinterval  of  the  total  range  of  the  independent  variable.  The  square  of 
the  ratio 


y  =  |H-HR|/S  (40) 

appears  in  (9)  due  to  the  fact  that  we  have  selected  the  function  to  be  (1-x)*/2 
where,  of  course,  x  is  equivalent  to  y2.  We  point  out  that  if  one  chooses  the 
function  to  be  (1-y2)^/2  an<j  proceeds  with  a  development  similar  to  that  of  the 
preceding  paragraphs,  then  the  result  will  be  an  algorithm  along  the  lines  of  (7)  - 
(9)  with  the  exception  that  |H-Hr|/S  replaces  its  square  in  (9)  and  yn  =  J  in  (7). 
However,  elimination  of  the  squaring  operation  exacts  a  price  in  terms  of  memory. 
In  particular,  the  curvature  of  the  graph  of  the  function  (1-y2)^'2  is  greater 
than  that  of  (1-x)^2  on  the  interval  from  0  to  1.  Consequently,  for  a  given 
error  bound  e,  one  can  expect  that  a  greater  number  of  straight  line  segments  will 
be  required  to  approximate  the  former  function.  For  instance,  if  is  1/8  nmi, 
then,  as  shown  in  figure  2,  the  algorithm  (7)  -  (9)  requires  a  partition  of  the 
interval  from  0  to  J2  consisting  of  four  subintervals.  On  the  other  hand,  if  the 
squaring  operation  is  eliminated,  then  it  can  be  shown  that  this  error  bound  can  be 
satisfied  by  a  similar  subdivision  of  the  interval  from  0  to  J  provided  that  the 
number  of  subintervals  is  at  least  16. 


7.  CONCLUDING  REMARKS 


We  have  introduced  a  systematic  approach  to  the  design  of  algorithms  for  estimating 
corrected  slant  range,  subject  to  prescribed  limitations  on  the  estimation  error. 
In  terms  of  accuracy,  utilization  of  memory,  and  computational  speed,  the  approach 
is  capable  of  providing  an  algorithm  that  is  superior  to  the  method  for  determining 
corrected  slant  range  that  is  currently  employed  in  NAS.  It  also  provides  a 
vehicle  for  meeting  more  stringent  error  requirements  in  the  future  that  are  con¬ 
sistent  with  the  needs  of  advanced  air  traffic  control  functions  such  as  Conflict 
Alert . 
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APPENDIX  A 


PROPERTIES  OF  THE  FUNCTION  E(u,0) 

Suppose  u  is  a  member  of  the  continuum  from  0  to  J  2 ,  and  let  us  consider  E(u,0) 
as  a  function  of  0  alone.  From  (18)  -  (21)  and  (24),  it  is  clear  that  the  function 
vanishes  at  0=  u.  Other  characteristics  of  the  function  can  be  determined  by 
examining  the  partial  derivative  with  respect  to  0.  First,  it  can  be  shown  that 
the  function  monotonical ly  decreases  as  0  increases  from  0  to  u,  and  thereafter  it 
monotonic al ly  increases  as  0  moves  from  u  toward  J2.  Moreover,  if  uj^u2,  then 
E(uj,0)  is  less  than  E(u2,8)  when  0<up  and  the  reverse  is  true  when  g  >u2-  As  a 
result,  the  functions  Efu^.B)  and  E(u2,8)  intersect  at  one  and  only  one  point  in 
the  continuum  from  0  to  J2  and  this  is  a  number  0  that  is  greater  than  u^  and  less 
than  u2-  The  situation  is  illustrated  in  figure  A-l . 

We  now  point  out  that  when  the  function  E(u,  0)  is  viewed  as  a  function  of  the 
vector  (u,0)  on  the  set  of  values  for  which 

0<u<J2  and  O<0<J2  (A-l) 

it  is  continuous.  Moreover,  because  the  set  is  closed,  it  is  uniformly  continuous 
(see  reference  3).i.e.,  to  each  e>0  there  corresponds  some  e >0  such  that 

jEtvpWj)  -  E  (v2,w2)|_<e  (A-2) 

whenever  the  Euclidean  distance  between  (vj.Wj)  and  (v2,w2  )  is  less  than  a. 
Thus,  recognizing  that  E(u,0)  vanishes  at  0=u,  it  follows  directly  from 

figure  A-l  that  to  any  error-bound  e>0  there  corresponds  a  positive  integer  mt  such 
that 

E(ult0*)  =  E(u2.8*)  (A-3) 

whenever  m  >  m  and 

-  O 

0  iu2  '  ul<2~m  (A-4) 

Suppose  now  that  u  lies  between  u^  and  u2  of  figure  A-l.  Then,  from  the  preceding 
paragraph,  we  know  that 

/  E(uj , 0*)  if  u  <  0* 

E(u,0*)<  /  (A-5) 

^  E(u2, 0*)  if  u  >  0* 

Consequently,  E(u,  0*)  cannot  exceed  c  whenever  m>m4  and  uj  and  u2  satisfy (A-4). 
The  statements  made  in  connection  with  (29)  and  (30T  follow  directly. 
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APPENDIX  B 


ERROR  BOUND  AFTER  p  NEWTON-RAPHSON  ITERATION 


Using  (38),  it  can  be  shown  by  direct  calculation  that 

a 

zk+l  *  *  "  <2k  _a)2  /2* k  (B_1) 

Moreover,  since  the  first  and  second  derivatives  of  the  function  (37)  with  respect 
to  z  are  positive  for  all  z>0,  it  follows  that  the  tangent  to  the  graph  of  the 
function  at  any  point  z^O  mu8C  intersect  the  horizontal  or  z-axis  at  a  point  z^+i 
greater  than  or  equal  to  a •  In  other  words,  starting  with  an  initial  estimate 
z  >0  of  corrected  slant  range,  the  estimates  zpZ2,...  provided  by  the  Newton- 
Raphson  iteration  will  all  be  at  least  as  large  as  a-  Consequently,  if  a  is  known 
to  be  greater  than  1/2  nmi,  2  za  is  greater  than  1  mi,  and  [  z(  -  «  |  does  not 
exceed  e/ 2  nmi,  then  (39)  follows  directly  from  (B-l)  for  all  integers  p >1 .  As 
will  be  shown  next,  all  three  of  these  conditions  can  be  met  under  typical  NAS 
operations . 

Relations  (13)  and  (15)  -  (17)  imply 

«  -  R(S,x)  >Sm  (1-J2)1/2  (B-2) 


Also ,  since  Ea  (x)  does  not  exceed c  for  the  estimate  (32)  provided  by  the  algorithm 
(7)  -  (9),  it  follows  from  (34)  and  (36)  that  the  absolute  value  of  the  difference 
between  zq  and  a  cannot  be  greater  than  c/2.  As  a  result,  z#  will  never  be  less 
than  «-  c/2.  Thus,  from  (B-2),  we  conclude  that 

2z#  >  2a -c^  2Sm  (l-j2)1/2  -c  (B-3) 

But  the  right  side  of  (B— 3 )  is  greater  than  or  equal  to  1  if 

e  <2  Sm  (l-J2)l/2-i  (b-4) 

Hence,  when  Sm  is  2  nmi  and  J  is  sin  70*,  (B-2)  implies  cr  is  greater  than  1/2,  and, 
from  (B-4),  it  follows  that  2z  must  exceed  1  if  the  error  bound  c  is  less  than  4 
cos  70*  -1  (0.368). 
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